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■i-H . 1 INTRODUCTION 



ABSTRACT 

The Karhunen-Loeve (KL) eigenvectors and eigenvalues of the sample correlation ma- 
trix are used to analyse the spatial fluctuations of the REFLEX clusters of galaxies. 
The method avoids the disturbing effects of correlated power spectral densities which 
affects all previous cluster measurements on Gpc scales. Comprehensive tests use a 
large set of independent REFLEX-like mock cluster samples extracted from the Hubble 
Volume Simulation. It is found that unbiased measurements on Gpc scales are possible 
with the REFLEX data. The distribution of the KL eigenvalues are consistent with a 
Gaussian random field on the 93.4% confidence level. Assuming spatially flat cold dark 
matter models, the marginalization of the likelihood contours over different sample 
volumes, fiducial cosmologies, mass/X-ray luminosity relations and baryon densities, 
yields the 95.4% confidence interval for the matter density of 0.03 < ft m h 2 < 0.19. 
The N-body simulations show that cosmic variance, although difficult to estimate, is 
expected to increase the confidence intervals by about 50%. 
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The cosmological parameters characterize the time evolution 
of the cosmic scale factor, and determine the formation and 
evolution of structures within the Universe. Rich clusters of 
galaxies are physically well-defined tracers of these struc- 
tures because they can only be formed at well-defined sites, 
namely where the peaks of the initial density field exceed 
a critical density threshold. This threshold is soley deter- 
mined by gravitation. Gaussian initial conditions simplify 
the situation even more. Therefore, the physical properties 
of the cluster population, like mass function and spatial dis- 
tribution, are closely related to the global properties of the 
Universe and give thus direct information on the values of 
the cosmological parameters. 

Important constraints on the values of the cosmologi- 
cal parameters obtained with galaxy clusters are generally 
based on measurements of the mean cluster abundance (e.g., 
Viana & Liddle 1996, Bahcall & Fan 1998, Borgani et al. 
2001, Reiprich & Bohringer 2002, see also the theoretical 
work of Haiman, Mohr & Holder 2001). However, the clus- 
ter abundance probes only a small scale range so that the 
resulting values of the matter density and the normalization 



parameter, erg, of the structure formation models are highly 
correlated. 

Measurements of the spatial fluctuations of the cluster 
abundance over a sufficiently large scale range can break the 
degeneracy. A review of recent measurements obtained with 
the spatial two-point correlation function of galaxy clusters 
is given in Collins et al. (2000). The fluctuations are also 
characterised by the power spectrum, P(k), which is di- 
rectly related to theory. Recent measurements of this quan- 
tity use either optically selected clusters (Peacock & West 
1992, Einasto et al. 1993, Jing & Valdarnini 1993, Einasto 
et al. 1998, Retzlaff et al. 1998, Tadros et al. 1998, Miller 
& Batuski 2001) or X-ray selected clusters (Retzlaff 1999, 
Schuecker et al. 2001, Zandivarez, Abadi & Lambas 2001). 
The advantages of X-ray against optically selected cluster 
samples are discussed in, e.g., Borgani & Guzzo (2001). 

For the construction of the ROSAT ESO Flux Limited 
X-Ray (REFLEX) cluster sample special care was taken 
to get a homogeneous sampling and a high completeness 
(Bohringer et al. 2001). The sample consists of 452 clus- 
ters with redshifts z < 0.45, selected in X-rays from the 
ROSAT All-Sky Survey and is confirmed by extensive op- 
tical follow-up observations within a large ESO Key Pro- 
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gramme (Bohringer et al. 1998, Guzzo et al. 1999). This 
makes the sample well-suited for spatial analyses on Gpc 
scales. 

However, on Gpc scales the anisotropy of the volumes 
of all cluster surveys becomes apparent. Therefore, reliable 
P(k) measurements of the projects mentioned above could 
only be obtained up to maximum scales reaching 200 to 
400 h~ x Mpc. Unfortunately, the plane waves used in the 
standard power spectrum analyses to expand the observed 
fluctuations are no longer orthogonal on Gpc scales and must 
be replaced by another set of basis functions fulfilling this 
fundamental criterion. The conditions of orthogonality of 
the basis functions and statistical orthogonality of the ex- 
pansion coefficients lead to the Karhunen-Loeve (KL) eigen- 
vectors of the sample correlation matrix (Karhunen 1947, 
Loeve 1948). They offer an analysis of the cluster power 
spectrum which is free from any disturbing effects of corre- 
lated power spectral densities affecting all previous cluster 
measurements on Gpc scales. 

The present paper applies the KL method to estimate 
the cosmic matter density and the linear normalization of 
the matter power spectrum using the spatial fluctuations of 
the REFLEX clusters. In order to introduce the basic quan- 
tities and to make the paper more self-contained, we recall 
in Sect.^ some aspects of the KL method and its applica- 
tion to large-scale structure work. The relations between the 
observed quantities as measured in the present investigation 
and the cosmological parameters are derived in Sect.fel The 
basic properties of the REFLEX cluster sample are sum- 
marized in Sect.^j. The KL eigenvectors and the spectrum 
of the eigenvalues of the REFLEX sample are presented in 
Sect.^. The final results on the cosmic matter density and 
normalization of the matter power spectrum obtained with 
the REFLEX sample are given in Sect.|^ and are discussed 
in Sect.f?]. 

To evaluate systematic and statistical errors as well as 
the effects of cosmic variance, end-to-end tests are performed 
which follow the basic steps of the REFLEX survey reduc- 
tion and the KL method of parameter estimation. Here we 
use a large set of independent REFLEX-like mock cluster 
samples selected from the Hubble Volume Simulation. The 
details are given in Appendix |a| 

As the fiducial cosmological model which is used to com- 
pute geometric quantities and KL eigenvectors, we assume a 
pressure-less, spatially flat Friedmann-Lemaftre model, the 
cosmic matter density, fi m = 0.3, the cosmological constant 
in the form Qa = 0.7, and the Hubble constant in units of 
h = i/o/lOOkms -1 Mpc -1 . 



the power spectrum (cosmological) parameters which maxi- 
mizes the probability of the observed fluctuations (Sect. 2.2) 



2.1 Calculation of the eigenvectors 

The survey volume is devided into M cells, each with a spe- 
cific comoving volume, Vi. We chose spherical coordinates 
and specified each edge of a cell by the three normal Euler 
coordinates. The results of the KL analysis do, however, not 
depend on a specific pixellation (see below). 

In the i-th cell centered on the comoving coordinate 
vector ft, Di clusters are counted. The expansion of the 
field, Di, can be written in the component form as Di = 
Y^jLi Bj , i = 1, . . . , M, where the tpij are the ele- 
ments of a matrix which gives the i-th component of the 
j-th basis vector. 

The modes and coefficients should fulfill two crite- 
ria, (i) The modes should be orthogonal to each other, 
4>Tk' l Pkj = Sij, where T denotes the transpose of a 
matrix and Sij the Kronecker delta, (ii) The modes should 
yield statistically orthogonal expansion coefficients. One 
thus requires that the expectation value of the sample co- 
variance matrix has the form < Bi Bj > = < Bj > Sij . The 
two criteria directly lead to the equations which determine 
the optimal basis vectors, 



^ Rki ipji = < Bj > Tpkj = Xj i>kj ■ 



(1) 



with the components of the correlation matrix, Rki, 
defined via the expectation values, < Bi Bj > = 

J2t*,l=i i'JkRkli'lj, through Ru =< D k Dj >. The prob- 
lem of finding the set of modes satisfying the conditions 
of orthogonality and statistical orthogonality thus reduces 
to the problem of finding the eigenvectors of the correlation 
matrix R, called the KL eigenvectors, and the corresponding 
eigenvalues, Aj =< Bf >, constituting the KL fluctuation 
spectrum. 

For an arbitrary pixellation of the survey volume 
the noise per counting cell varies even for volume-limited 
samples, and one has to diagonalize the noise compo- 
nent of the correlation matrix before the eigenvectors 
are computed. The separation of signal and noise in 
the new basis is achieved by transforming (whitening) 
the elements of the correlation matrix computing R^j — 

Y^k,i=i N ik 1/2 R kl N ij 1/2 - Tne N u: 1/2 are tne inverse square 
roots of the elements of the noise correlation matrix, Nik = 
8ik Jy < n(r) > d 3 r = Ni, and < n(f) > the expected 
cluster number density at the comoving position r. 



2 THE KL METHOD 

The KL method was first used to test cosmological structure 
formation models by Bond (1995) using cosmic microwave 
background (CMB) temperature maps. Vogeley & Szalay 
(1996) translated the method to the case of the spatial anal- 
ysis of galaxy distributions. Applications to galaxy surveys 
can be found in Matsubara, Szalay & Landy (2000) and Sza- 
lay et al. (2001). The KL method as used here to analyse 
cluster data cons ists of two steps: calculation of the eigen- 
vectors (Sect. 2.1), and likelihood estimation of the values of 



2.2 Estimation of model parameters of the power 
spectrum 

The present investigation tests the fluctuating part of the 
cluster number counts. Therefore, the covariance matrix, 
C, of the KL coefficients is used to estimate the values of 
the (cosmological) parameters, x\, . . . , x q , characterizing the 
power spectrum. 

The covariance matrix is estimated in the following way. 
Choose a specific set of Xi values to specify the model P(k). 
Fourier-transform P(k) by direct numerical integration in 
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order to get the correlation function, £(r), and compute the 
continuous part of the cell-averaged correlation matrix (we 
use a Monte-Carlo estimate) of the model to be tested: 



dVi / d A r 3 a\n-? 3 \). 



(2) 



Choose also an appropriate model for the expected aver- 
age number of clusters, N%, in each cell. To be consistent 
with the fluctuation analysis we use an empirical model (see 
Sect. 3.3). This model is not changed during the testing of 



different model power spectra. The coefficients 



i 



Al 



(N k Ni$ k , + N k Su) 



jhj_ 



(3) 



constitute the estimated model covariance matrix, C, of the 
KL coefficients, where the first term on the right-hand side 
of (^) describes the clustering signal and the second term 
the noise. The ipij are the KL eigenvectors of the whitened 
corre latio n matrix obtained with the fiducial cosmology (see 
Sect. 2.1). The model covariance matrix is not diagonal un- 
less the fiducial model used to compute the KL eigenvectors 
is identical to the model used to compute the 

In Sect.[] it will be shown that the frequency distribu- 
tion of the REFLEX KL coefficients, Bi, is well described 
by a Gaussian. Due to the linearity of the KL transform this 
suggests that the REFLEX cluster density field is governed 
by a Gaussian-like random field (for large cell sizes). The 
multivariate likelihood function of the parameters Xi should 
thus be of the form 



C(Bi, . . . , Bm\xi, . . . ,x q 
(27r)- M/2 |detCr 1/2 



. x]» I -^AB T C' 1 AB 



(4) 



with the difference vector AB = B— < B >. The Xi val- 
ues of the power spectrum parameters which maximise the 
probability of obtaining fluctuations transformed into the 
KL base as large as observed are defined by the maximum 
of the sample function (0). 



3 THE RELATIONS BETWEEN THE 
OBSERVED QUANTITIES AND THE 
COSMOLOGICAL PARAMETERS 

The KL method estimates the values of the cosmological pa- 
rameters comparing the observed fluctuations of the cluster 
number densities transformed into the KL eigenvector basis 
with theoretical expectations. In the following the ge nera l 
assumptions on the matter power spectrum (see Sect. 3.1), 
on the relation between the observed amplitude of the clus- 
ter power spectrum and the standar d as normalization of 
the matter power spectrum (see Sect. 3.2), and on the em- 
pirical model used to compute the expected mean cluster 
number counts, Ni, (see Sect. 3.3) are described. 



3.1 Matter power spectrum 

On the largest scales the density field is assumed to be Gaus- 
sian with a power-law spectrum of adiabatic matter fluctu- 
ations, P(k) ~ k n . In order to describe the matter power 
spectrum on smaller scales, we take into account the effects 



of a collisionless matter component and the collisional bary- 
onic component. Instead of solving the corresponding mul- 
tispecies Boltzmann equations for each model to be tested, 
the comparatively simple fitting formulae for the transfer 
functions, T x (k), given in Eisenstein & Hu (1998) are used 
providing a more accurate description than the standard 
BBKS fitting functions mainly characterized by the scale- 
independent shape parameter, F. 

We restrict the present KL analyses to the estimation 
of the matter density, Q m = x and as because they de- 
termine - for a given Hubble constant - the general shape 
and amplitude of the power spectrum. For the Hubble con- 
stant we take h = 0.7 as suggested by Hubble Space Tele- 
scope observations (Freedman et al. 2001). The final re- 
sults on £l m are given in units of h. In addition, we assume 
n — 1.0, a mean temperature of the CMB of Tcmb = 2.728, 
a spatially flat cosmology as suggested by CMB measure- 
ments (De Bernardis et al. 2000), and the baryon density 
Q b h 2 = 0.0196 as suggested by chemical abundance mea- 
surements of distant quasars (Buries & Tytler 1998) and 
Standard Big Bang Nucleosynthesis calculations (Buries, 
Nollett & Turner 2001). In the next paper additional ob- 
servational constraints will be included so that the priors 
can be weakened. 



3.2 Relation between as and the observed 
amplitude, Po, of the power spectrum 

It is generally assumed that on large scales structure growth 
can be treated within linear theory. The observed cluster 
power spectrum, P b s (k), is the result of a complex averag- 
ing process over evolving matter power spectra, P(k, z) = 
P(k)D 2 (z) I D 2 (0) and clusters with different values of the 
biasing parameter, b(M,z). Here, D(z) is the linear growth 
factor. Matarrese et al. (1997) and Moscardini et al. (2000) 
derived analytic equations approximating this process which 
we summarize by the equation 



Pobs(fc) = (P(k,z) (b(M,z)) 



m) z 



(•>) 



where the mass and redshift expectations involve the actual 
number of clusters, N(M, z)dMdz, observed within given 
mass and redshift shells, and the corresponding redshift his- 
togram, N(z)dz. Here, < b(M, z) >m is the mass-weighted 
biasing factor and < ■ >z the redshift average (eq. 14 in 
Moscardini et al. 2000). Within the general framework of 
linear perturbation theory of cosmic structures, the present- 
day matter power spectrum is 



P(k) 



2tv 2 aj k n T 2 (k) 



J dfcfc"+ 2 T|(fc)|VK(8fc//iMpc- 1 )| 2 ' 



(6) 



P(k) = a\ k n T*(k) Cnig, (7) 
where we have introduced for convenience the function 

(8) 



(nxR = ^2 / T 2 (k)\W(kR)\ 2 



The spectrum is normalized by as in the standard way using 
the Fourier-transformed top-hat filter, W(kR), with the co- 
moving radius R = 8 ft -1 Mpc. It is important to note that 
as defined in this way reflects the amplitude of the power 
spectrum without any non-linear corrections. 
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Figure 1. Examples of one-dimensional tracings of the three-dimensional KL eigenvectors as a function of Right Ascension (left), 
Declination (middle), and comoving radial distance (right). The values of the eigenvectors are computed for each direction (in arbitrary 
physical units) at the centers the 10 cells and are then interpolated (for illustration) by cubic splines. Continuous lines show eigenvectors 
with the highest eigenvalue, dashed lines with lower eigenvalues (higher orders). The REFLEX survey crosses the galactic plane at 
R.A. = 120 deg and 270 deg, and has no clusters in the North as seen by the low values of the eigenvectors at the corresponding positions. 
The effective depth of the REFLEX survey is at r = 150 h" 1 Mpc. 



Equation (g) implies that the observed mass and 
redshift-averaged cluster power spectrum has the same 
shape, i.e., functional form as the underlying matter power 
spectrum, P(k). Therefore, we set P ohB (k) = P k n T 2 (k) 
with the parameters x and Po determined by observation. 
Equating the latter formula and the theoretical expectation 
(^), yields the general relation between the observed ampli- 
tude and the normalization of the matter power spectrum, 



2 f-\ 
a S Smsct 



D 2 (0) 



(b(M,z)) 



(9) 



As expected, the observed amplitude, Po, of the power spec- 
trum depends on the sample. Note that in the present case, 
the actual values of the cluster sample are inserted in (||), 
where the masses are obtained from the observed X-ray lu- 
minosities using the empirical mass-to-X-ray luminosity re- 
lation (see eq. |Al|) of Reiprich & Bohringer (2002). 

Equation (H) shows that a specific biasing model has to 
be chosen in order to deduce from Po the linear normaliza- 
tion, as- 



3.2.1 High-peak biasing 

The KL likelihood analysis is based on the assumption of 
a Gaussian random field, supported by the observed distri- 
bution of the REFLEX KL eigenvalues (see Sect.^). In the 
line of this observation we apply the related biasing scheme, 
b(M,z) = 8 c (z)/a 2 (M,z), derived by Kaiser (1984) for 
galaxy clusters on the same statistical grounds as a Gaussian 
random field in the limit of high density peaks. The two con- 
ditions, a{M,z) « 8 c (z), and (5 c (z)/a 2 (M, z)) 2 £(r, z) < 1, 
are generally fulfilled in the present case because the fluc- 
tuation analyses are performed with massive clusters where 
l<&<4on50<r< 1000 ft -1 Mpc scales where the mat- 
ter correlation function is about 10~ 2 > |£| > 10 -6 . The 
results obtained with the simulations shown in Appendix^ 
are also consistent with these assumptions. In the Kaiser 
model the biasing parameter is determined by the slightly 
redshift-dependent critical density threshold, 5 c (z), of the 



spherical collapse model, and the mass variance, a 2 (M, z). 
The latter quantity can be written in terms of £ as 



a\M{R),z) = *i 



2 D 2 (z) Cn 



(10) 



D 2 (0) ' 

The variance in eq. ([lo]) decreases with z in a manner that 
at high redshift the biasing for clusters of a given mass is 
stronger than the decreasing matter power spectrum. For 
galaxy clusters, Po is thus expected to increase with z. In- 
dependent from any redshift-dependent effect, the high-peak 
biasing gives the monotonic relation Po ~ l/crf- 



3.2.2 Other biasing schemes 

Based on the Press-Schechter prescription, Mo & White 
(1996) derived a formula which describes the biasing of 
galaxy-size objects, 



6(M, z) 



5 c {z) 1_ 

a 2 (M,z) S c (z) 



(11) 



As for the high-peak biasing, the Mo & White biasing de- 
pends via the second term on the right-hand side of ( |ll| ) on 
as, however, now with two additional terms. The first term 
describes the peculiar motions of the dark matter haloes 
and the second and third terms the effects of the peak- 
background split. In contrast to the high-peak biasing, the 
relation between <rf and Po has a quadratic form. Therefore, 
the model suggests two values of erg for a given Po. The high 
as case characterizes an almost unbiased halo distribution 
where basically each mass peak corresponds to a virialized 
object (low-biasing regime, unrealistic case for clusters). The 
small as case characterizes a strongly biased distribution 
where the virialized structures appear as rare objects (high 
biasing regime, realistic case for clusters). 

The biasing formula given in Sheth & Tormen (1999) 
has the same properties as (|ll]). It is, however, better cali- 
brated with N-body simulations over a mass range reaching 
5 10 13 h' 1 M Q (or the X-ray luminosity 6 10 42 h~ 2 erg s" 1 for 
the energy range 0.1 — 2.4 keV using eq. Al). Unfortunately, 
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counts required for a safe detection, the local satellite's ex- 
posure time, and the local galactic neutral hydrogen column 
density, iVm (see Sect.^ and Bohringer et al. 2001). Random 
number distributions are computed to generate angular dis- 
tributions which precisely follow the survey boundaries as 
described in Collins et al. (2000) and Schuecker et al. (2001). 

For the radial part of Ni we also generate random dis- 
tributions which are now guided by the observed redshift 
histogram smoothed with a Gaussian filter profile with the 
standard deviation cr z = 0.015. We compared the KL like- 
lihood contours obtained with the smoothing method and 
with the X-ray luminosity function given in Bohringer et al. 
(2002). No significant differences are found in the values of 
the estimated parameters as long as a significant number of 
clusters have comoving distances reaching > 300 h" 1 Mpc. 



Figure 2. The spectrum of the KL eigenvalues, < B„ >, as the 
function of rank of the REFLEX cluster sample obtained for the 
fiducial cosmology. 



150 



100 




-2024 
(B n -<B n >)/V<B=> 

Figure 3. Histogram of the normalized KL eigenvalues (mean 
0.059, standard deviation 0.968) with superposed normal Gaus- 
sian profile for the REFLEX cluster sample. 



this maximum X-ray luminosity is close to the minimum X- 
ray luminosity for completeness of the cluster sample used 
in the present investigation (see Sect. HI). 



3.3 Empirical model for the average cluster 
number densities 

The present investigation concentrates on the exact mod- 
elling of the fluctuating part of the cluster number counts. 
The model for the average cluster number densities, Ni (as 
used in eq. |j) , is thus not changed during the likelihood op- 
timization. The Ni are estimated by the following empirical 
Monte-Carlo method. 

For an X-ray cluster sample, the angular part of Ni is 
mainly determined by the local X-ray flux limit of the sur- 
vey, which in turn is given by the preset nominal flux limit 
of the sample (see Sect.M), the minimum number of source 



4 THE REFLEX SAMPLE 

The REFLEX sample has 452 southern X-ray clusters of 
galaxies, 449 with measured redshifts, z < 0.45 (Bohringer 
et al. 2001). The clusters are selected in an area of 13 924 
square degrees (4.24 sr) from the ROSAT All-Sky Survey 
(RASS, Triimper 1993, Voges et al. 1999). The nominal limit 
of the unabsorbed X-ray fluxes is 3 10~ 12 erg s -1 cm" 2 in the 
energy range (0.1 — 2.4) keV. 65 percent of the sample are 
Abell/ACO/Supplement clusters. 

In order to reduce strong spatial variations of the sam- 
pling, the 452 REFLEX clusters were selected outside the 
galactic plane (galactic latitudes |fo| > 20deg) and some ad- 
ditional crowded stellar fields (e.g., Magellanic Clouds). The 
remaining corrections for the satellite exposure time and 
galactic absorption are well-documented and can be mod- 
elled in detail (e.g., Bohringer et al. 2001). The sample has 
been successfully used for the determination of the X-ray 
luminosity function (Bohringer et al. 2002), for the analyses 
of the cluster correlation function (Collins et al. 2000), the 
related peculiar motions (L. Guzzo et al., in preparation), 
and the power spectrum (Schuecker et al. 2001). 

Several incompleteness tests described in these RE- 
FLEX papers are based on either the REFLEX sample it- 
self or other observed or simulated cluster samples. The 
tests suggest the absence of a significant incompleteness for 
clusters with X-ray luminosities Lx > 2.5 10 42 /i -2 ergs -1 . 
The present investigation uses the 428 clusters which have 
at least 10 X-ray source counts and which fall within this 
well-controlled luminosity range. For comoving distances 
r < 500 /i -1 Mpc (z < 0.18) no systematic deficiencies of 
the comoving REFLEX cluster number densities are found. 



5 THE REFLEX KL EIGENVECTORS AND 
EIGENVALUES 

A spherical volume containing the REFLEX survey up to a 
certain maximum comoving radius, r, is devided into 1000 
volume elements (spherical coordinates): 10 angular bins in 
Right Ascension, 10 in Declination, and 10 bins along the 
comoving radial axis. The number s of REFLEX clusters and 
random sample points (see Sect. 3.3) in each of the cells 
are counted, and standard linear algebra codes (Press et 
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Figure 4. Likelihood contours (68.3, 95.4, 99.0%) of the REFLEX sample in the Po-Q m (left panels) and in the ag-Q m parameter space 
(right panels). The erg values are computed with the high-peak biasing (Kaiser 1984). The amplitudes, Pq, of the power spectra are given 
in units of h -4 Mpc 4 . The upper row show the results for the 342 REFLEX clusters located within comoving distances r < 500 Mpc -1 
(z < 0.18) for the fiducial cosmology. The lower row shows the results for the 403 clusters within r < 750 h~ Mpc (z < 0.27) for the 
same fiducial cosmology. The crosses mark the points with the highest likelihood value. 



al. 1989) are used to compute the eigenvectors and eigenval- 
ues of the pixel-averaged whitened correlation matrix. The 
following KL analysis is restricted to the M = 540 eigenvec- 
tors with nonzero eigenvalues, sorted (ranked) with decreas- 
ing < B% >, i.e., with decreasing signal-to-noise. A few ex- 
amples of one-dimensional tracings of the three-dimensional 
eigenvectors are shown for the largest eigenvalues in Fig.|l|. 

The spectrum of the REFLEX KL eigenvalues is shown 
in Fig. ^| The spectrum basically follows a power law. 
This indicates that, excluding the extreme n ranges, the 
KL eigenvectors sample three-dimensional structures over a 
large (but not the complete) n range. 

The frequency distribution of the normalized KL eigen- 
values gives information about the Gaussianity of the dis- 
crete fluctuation field and is thus quite important for the jus- 
tification of the multivariate Gaussian likelihood functions 
which will be used fo r the estimation of the power spectrum 
parameters (see Sect. ^21). The histogram of the normalized 
deviations shown in Fig.tfl is consistent with a Gaussian ran- 
dom distribution on the 93.4% confidence level (KS test) and 
thus supports our basic assumption. Note that this result is 
mainly determined with cells larger than (50 h~ x Mpc) 3 . For 
smaller cells deviations from Gaussianity are expected and 
other likelihood functions must be used. The linearity of the 
KL transform suggests that the Gaussian distribution of the 
KL coefficients translates into a Gaussian random field of 
the underlying matter distribution. This fa vours the biasing 
model proposed by Kaiser (1984, see Sect. 3.2.1 ). 



6 RESULTS 

The KL method was tested with 27 independent mock clus- 
ter samples selected from the Hubble Volume Simulation. 
The simulations have the same fiducial cosmology as used 
here, fi m = 0.3 and ag = 0.9, inclu ding the values of n, h, 
Qbh 2 mentioned in Sect. 



3.1 



The details are given in 



Tomb, 

Appendix |a| For each cluster sample and biasing model, f2„ 
and as are varied within suitable intervals, and the resulting 
model covariance matrixes Q are computed. The maximum 
of the likelihood (^) is used to select the best estimate of Q m 
and as for each sample and biasing model. The sample-to- 
sample variations of the parameter values give at least for 
the fiducial cosmology an estimate of the errors of Q m and 
as when the effects of cosmic variance are included. 

The mean and la errors as obtained from the simula- 
tions are for the matter density £l m = 0.28 ± 0.14 and for 
the linear matter normalization as = 0.87T0.32 (high-peak 
biasing), a s = 1.20T0.66 (Mo & White biasing), 0.82T0.43 
(Sheth & Tormen biasing). Note that for the computation 
of the mean and standard deviation of the latter two biasing 
models only the likelihood maximum w hich i s located in the 
high biasing regime was used (see Sect. 3.2.2). Below we will 



compare these errors with the errors provided by the like- 
lihood contours computed with the KL method. Compared 
to the input values of the Hubble Volume Simulation, no 
significant systematic errors are thus found (see also Fig. Al 
in Appendix |X|). 
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Figure 5. Upper row: Likelihood contours as in the upper row of Fig. [ij but for the Einstein-de Sitter fiducial cosmology. Lower row: 
Likelihood contours using the empirical mass/X-ray luminosity relation of Reiprich & Bohringer (2002, see eq. |Al| ), but with cluster 
masses artificially boosted by a factor of two with respect to the observed relation. 



The KL method was then applied to the REFLEX clus- 
ter sample. The main results are plotted in the upper panels 
of Fig. ^| Shown are the 1-3<j likelihood contours in the Po- 
Qm and <7g-f2 m parameter spaces. The cosmic matter density 
with the highest likelihood value is Q m — 0.16 ± 0.06 (la 
error without cosmic variance and no marginalization with 
respect to h). It will be seen that the f2 m values are basically 
unaffected by the assumed biasing model used to compute 

The as values shown in the right panel of Fig.^J are 
determined with the biasing scheme of Kaiser (1984). The 
results obtained with the Mo & White (1996) and Sheth 
& Tormen (1999) biasing models are discussed below (see 
Fig. |^) . The linear normalization of the matter power spec- 
trum with the highest likelihood value is as = 1.2 ± 0.3 (la 
error without cosmic variance and no marginalization with 
respect to h). 

The sensitivity of the results on several given parameter 
values is illustrated in the lower panels of Fig. ^| and in Figs. ^ 
to[| 

In Fig. ^ the results obtained within a maximum comov- 
ing distance of r = 500 hT 1 Mpc corresponding to z = 0.18 
(upper panels) are compared to the results obtained within 
r = 750 h" 1 Mpc or z = 0.27 (lower panels). It is seen that 
not much information on the spatial fluctuations is gained 
by the KL method when the REFLEX clusters outside the 
well-tested redshift range (see Sect.^J) are included. 

The upper panels of Fig. ^| show the likelihood contours 
determined with an Einstein-de Sitter fiducial cosmology. 
Due to the fact that the large-scale structures are mainly 



probed at z < 0.18, the effects of different fiducial cosmolo- 
gies are not very large and do not really modify the present 
KL results. 

In the lower panels of Fig. [| the KL results are shown 
where we used the empirical mass/X-ray luminosity relation 
of Reiprich & Bohringer (2002), but with a systematic shift 
applied towards larger X-ray masses by a factor of two, or 
equivalently a shift towards smaller X-ray luminosities by a 
factor of 2.5. The shift should "compensate" for several pos- 
sible sources of systematic errors which could modify the 
empirical relation (e.g., underestimated X-ray masses, con- 
tamination of the X-ray flux by active galactic nuclei, re- 
lations derived from flux-limited samples). The KL results 
show that even large changes in the mass/luminosity rela- 
tion in the given directions do not affect the estimation of 
the matter density. In the present case, only the normaliza- 
tion of the matter power spectrum is increased by 20%. 

In the upper panels of Fig. ^ we show the results ob- 
tained with the fiducial cosmology and the 2a upper limit 
on the baryon density of Qbh 2 = 0.029 obtained from the 
combination of BOOMERanG and COBE/DMR data (Masi 
et al. 2002). The main effect of the baryons is to steepen 
P(k). A large ttb value can thus be compensated by a large 
Q, m as seen in Fig.^. The same effect is found in the 2dF 
100k data (Percival et al. 2001, Tegmark, Hamilton & Xu 
2001). 

In Fig.^j we show the likelihood contour s de termined 
with all three biasing models described in Sect. 3.2. As men- 
tioned above, Q m is mainly independent of the assumed bi- 
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Figure 6. Likelihood contours for the fiducial cosmology but with the 2<r upper limit of the baryon density, Sl^h = 0.029, measured 
with the BOOMERanG experiment (Masi et al. 2002) plotted in the Po-Q m parameter space (upper left) and for the biasing models of 
Kaiser (1984, upper right), Mo & White (1996, lower left) and Sheth & Tormen (1999, lower right). Note the different scalings of the <rg 
axes. 



asing model, but an effect is seen in the derived erg values 
which will be discussed below in more detail. 

The most important cosmological constraint derived 
from the spatial fluctuations of the REFLEX clusters is the 
cosmic matter density obtained from the marginalization of 
the likelihood distributions shown in Figs.^]to|]. For h = 0.7 
the REFLEX data give the 95.4% confidence interval 



0.07 < ttm < 0.38 (95.4% without 



cosmic variance 



(12) 



Note that for fl(,/i 2 = 0.029 the highest likelihood value of 
most models is at Q m — 0.20. A more systematic analysis 
of models with different f2(, values (and n) is necessary and 
will be given in the next paper. 

The KL analysis of the spatial fluctuations of the RE- 
FLEX clusters is less sensitive to the linear erg. From the 
marginalization of the likelihood distributions and for the 
high-peak biasing model of Kaiser (1984) we obtained 



0.6 < ag < 2.6 (95.4% without cosmic variance) 



(13) 



with the highest likelihood value at as ~ 1.2. 

For the biasing schemes of Mo & White (1996) and 
Sheth & Tormen (1999) the situation is more complex. In the 
error-free case one would expect for the two biasing models 
two well-separated likelihood regions cente red o n the same 
fl m but at two different erg values (see Sect. 3.2.2). However, 



the comparatively large statistical scatter of the observed 
Po values smeares out the high- and low-biasing regimes and 
thus leads to the 'shoe-like' contours seen in the lower panels 
of Fig.|^. Nevertheless, the results obtained with the simula- 



Tab. 1. Comparison of the 95.4% confidence ranges for f2 m /i 2 ob- 
tained with galaxy clusters (REFLEX), recent measurements of 
CBM temperature fluctuations (BOOMERanG, DASI) and with 
galaxies (SDSS, 2dFGRS). References: (1) this work, (2) Netter- 
field et al. (2001), (3) Szalay et al. (2001), (4) Pryke et al. (2001), 
(5) Parcival et al. (2001). SDSS measures the shape parameter, 
r, 2dFGRS measures Q m h. These values are transformed using 
h = 0.7 and Q b h 2 = 0.0196. REFLEX, SDSS, and 2dFGRS as- 
sume a flat universe with a cosmological constant. BOOMERanG 
and DASI results have the weakest priors. 



Data 


Probe 


n 




Ref. 


REFLEX 


Clusters 


0.03 


- 0.19 


(1) 


BOOMERanG 


CMB 


0.05 


- 0.25 


(2) 


SDSS 


Galaxies 


0.08 


- 0.20 


(3) 


DASI 


CMB 


0.08 


- 0.24 


(4) 


2dFGRS 


Galaxies 


0.10 


- 0.18 


(5) 



tions (see Appendix ^) shows that all three biasing schemes 
give similar results, when for the Mo & White and Sheth 
& Tormen models the erg values located in the high-biasing 
regime are selected. 

The comparison of the errors of the parameter values 
obtained from the REFLEX data and from the simulations 
which include cosmic variance indicates that the errors in- 
cluding cosmic variance are about 50% larger compared to 
the KL errors. 
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7 DISCUSSION 

The present investigation applies the KL method to esti- 
mate the values of the cosmic matter density and the linear 
normalization of the matter power spectrum. The fluctua- 
tions of the comoving densities of the REFLEX clusters are 
analyzed up to Gpc scales with a well-defined survey specific 
set of eigenvectors. This offers the possibility to analyse the 
fluctuations up to Gpc scales without the disturbing effects 
of correlations between different power spectral densities, 
fobs(fc), which affects all previous cluster measurements on 
the largest scales. Note that the correlations artificially re- 
duce the statistical errors, so that simple numerical model 
fits to Pobs(fc) in order to estimate the values of the cosmo- 
logical parameters cannot be applied. 

The main result obtained with the KL analysis of the 
REFLEX clusters is that for spatially flat CDM-like struc- 
ture formation scenarios the data support a low-density uni- 
verse with (rescaling Q m to Q m h 2 using our prior h — 0.7) 

0.03 < Q m h 2 < 0.19, (14) 

and the linear normalization 0.6 < as < 2.6 (95.4% confi- 
dence intervals without cosmic variance). The notation un- 
derlines the fact that we did not marginalize the results for 
different h values. The Einstein-de Sitter case is ruled out 
with 99.99% confidence. The errors obtained with the KL 
method include marginalization over several important re- 
duction parameters but not cosmic variance. We have esti- 
mated the effect using 27 REFLEX-like mock samples se- 
lected from the Hubble Volume Simulation, and found that 
for the current sample the KL errors are probably underes- 
timated by 50%. 

We want to stress that the fi m measurements appear 
to be quite robust against several partially quite drastic 
changes of important reduction parameters. What really 
matters seems to be the baryon density (and thus also the 
spectral index, n, of the primordial power spectrum). A sys- 
tematic study of models with different and n is in prepa- 
ration. 

The REFLEX confidence range for Q. m h 2 in (0) is 
in good agreement with other recent measurements (see 
Tab. 1). The table gives the 95.4% confidence intervals ob- 
tained from different measurements with the minimum num- 
ber of priors (and not results obtained by combined data 
sets). Note that the different groups measured F, Q m , fl m h, 
or 0, m h 2 , and give either la or 2a errors. We have tried 
to transform the original results to Q m h 2 and 95.4% errors, 
having in mind that this can only be done approximately. 

The Sloan Digital Sky Survey (SDSS) result is obtained 
from the galaxy clustering of 222 square degrees early imag- 
ing data (Szalay et al. 2001). For SDSS the shape parameter, 
F, of the power spectrum is transformed to fl m h assum- 
ing h = 0.7 and Q,bh 2 = 0.0196 and the approximate for- 
mula given in Sugiyama (1995). The 2dF Galaxy Redshift 
Survey (2dFGRS) result for fl m h is obtained with 166 490 
galaxies. The Degree Angular Scale Interferometer (DASI) 
and BOOMERanG experiments measure the angular power 
spectrum of the CMB anisotropy. REFLEX, SDSS, and 
2dFGRS assume a flat universe with a cosmological con- 
stant, but the results do not strongly depend on S7a. They 
also assume n = 1. REFLEX has the additional constraint 
Q b h 2 = 0.0196. The BOOMERanG results have the weak 



prior 0.45 < h < 0.90 and eliminates models where the Uni- 
verse is younger than 10 Gyr. DASI assumes h > 0.45 and 
the optical depth due to reionization 0.0 < r c < 0.4. Com- 
pared to the results shown in Tab. 1, the REFLEX results 
extents to slightly smaller Qmh 2 values. Smaller confidence 
ranges from REFLEX are expected when the KL analysis 
will include both the fluctuations and the mean cluster num- 
ber densities, utilizing the complementarity of clustering and 
abundance of clusters. In this way the KL analysis will allow 
us to fully exploit the cosmological potential of the REFLEX 
survey of X-ray clusters. 

The 'banana-shape' likelihood contours obtained with 
the REFLEX data (see Figs. § to |) might be taken as an 
indication that the primordial power spectrum is less con- 
strained by the current REFLEX data. A significant im- 
provement is expected when the southern REFLEX sample 
and the northern NORAS sample (Bohringer et al. 2000) 
are extended to the deeper flux limit of 2 10 -12 erg s _1 cm -2 
in the energy range 0. 1-2.4 keV and combined to an all-sky 
sample of about 1 700 X-ray selected clusters of galaxies. 

We would like to thank the REFLEX group for their 
help in the preparation of the X-ray cluster sample, D. Eisen- 
stein and W. Hu for the computer code for the matter trans- 
fer functions, the Virgo Consortium for the simulated LCDM 
cluster sample, and the referee Stefano Borgani for his use- 
ful comments. P.S. acknowledges support under the grant 
No. 50 OR 9708 35. 
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APPENDIX A: VALIDATION OF THE KL 
ESTIMATION OF THE POWER SPECTRUM 
PARAMETERS WITH N-BODY SIMULATIONS 

Mock samples are used to test the likelihood method, es- 
pecially systematic errors and the effects of cosmic vari- 
ance. For studies of the clustering properties of X-ray se- 
lected cluster samples a crucial step is the transformation of 
the simulated cluster gravitational masses to the observable 
X-ray luminosities. Here we use the empirical mass/X-ray 
luminosity relation of Reiprich & Bohringer (2002) for the 
energy range (0.1-2.4) keV, 



/i" 2 10 44 ergs" 1 



= 7.199 x 10' 



M 



/) _1 Mq 



(Al) 



assuming a negligible intrinsic scatter. The formal la er- 
rors of the scaling factor and the index of (Al) are 6.3 and 



5.1%, respectively. In order to apply this relation, one has 
to ensure that the cluster masses as defined through the 
simulations are consistent with the masses of the empirical 
mass/X-ray luminosity relation as defined through the X- 
ray measurements. For the present error estimation we use 
the simple mass transformation model described below giv- 
ing redshift histograms similar to REFLEX. For exact model 



comparisons going beyond simple error estimates, more re- 
fined transformation models or simulations adapted to the 
empirical mass- luminosity relation should be used. 

Al Simulated clusters 

The Virgo Consortium provides cluster catalogues extracted 
from the Hubble Volume Simulations (see, e.g., Jenkins et 
al. 2001, see also Evrard et al. 2001). The public clus- 
ter catalogue used here is selected at z — from one 
ACold Dark Matter (LCDM) simulation with a box length 
of 3 7i _1 Gpc and fi m = 0.3, fi A = 0.7, and <j 8 = 0.90. 
The LCDM transfer function was computed with the CMB- 
FAST routine (Seljak & Zaldarriaga 1996) assuming h = 0.7, 
tt b h 2 = 0.0196 (Buries & Tytler 1998), and a primordial 
slope of P(k) of unity. Each of the 10 9 particles has a mass of 
2.22 10 12 h' 1 Mq. Motivated by the spherical collaps model 
the Virgo Consortium attempts to identify virialized regions 
that are overdense by a factor ~ 200 applying the friend- 
of-friend group finder with a linking length of 0.164. The 
resulting catalogue comprises 1 560 995 clusters. The mini- 
mum number of particles per cluster is 30. 

The friend-of-friend cluster masses obtained from the 



simulations are measured out to a radius. 



where the 



averaged density contrast relative to the local mass density 
is approximately 324. For the empirical mass/X-ray lumi- 
nosity relation, Reiprich & Bohringer used the radius r2oo, 
where the average density contrast of 200 is related to the 
Emstein-de Sitter critical mass density. The mass conver- 
sion factor is obtained from the relation between virial mass 
and X-ray temperature obtained from hydrodynamical sim- 
ulations (e.g., Bryan & Norman 1998), giving for fixed tem- 
perature and redshift M (r2oo) / M (r s i m ) = 0.69. A slightly 
better match between simulated and observed redshift his- 
tograms is obtained with the conversion factor 0.67 which 
we finally used. 

The M(r 200) cluster masses are transformed to X-ray 
luminosities in the energy range (0.1-2.4 keV). The observer 
restframe fluxes are obtained with the cluster luminosity 
distance taking into account the cosmic K-correction as ob- 
tained with a refined Raymond-Smith code (the cluster X- 
ray temperatures are estimated with the Lx-T relation of 
Markevitch (1998) without cooling flow corrections). The 
resulting total fluxes are reduced by 10 percent to get the 
measured fluxes because the X-ray observations do not in- 
clude the flux in the outer wings of the cluster X-ray image. 
This average difference between total and observed fluxes 
is obtained from Monte-Carlo simulations (H. Bohringer, in 
preparation). The variation of the X-ray flux limit of the 
REFLEX sample across the survey area are computed in 
the same way as in Schuecker et al. (2001). 

A2 Comparison of true and estimated parameter 
values 

We selected 27 independent REFLEX-like subsamples from 
the LCDM Hubble Volume cluster sample. The average 
number of clusters per sample and its standard deviation 
is 435 ± 28, similar to the 428 REFLEX clusters used for 
the final analyses. The redshift histograms closely resem- 
ble the observed distribution. We thus expect realistic error 
estimates from the simulations. 
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Figure Al. Histograms of the differences between the KL esti- 
mate of matter density, f2 m , and the linear normalization of the 
matter power spectrum, erg, with the input (true) values of the 
simulations. The biasing schemes are denoted by KA (high-peak 
biasing of Kaiser 1984), MW (Mo & White 1996) and ST (Sheth 
& Tormen 1998). The frequency distributions are obtained with 
27 REFLEX-like subsamples selected from the Hubble Volume 
Simulation. 



The histograms of the residuals between the values esti- 
mated with the KL method and the true (simulation input) 
values for Q m and for as as obtained w ith t he three biasing 
models (see Sect. 3.5) are shown in Fig. 



Al 



The frequency distribution of the KL estimates of (l m 
gives the formal mean and standard deviation of fi m = 
0.32 ± 0.19 (see Fig.0, upper left). Note that the distri- 
bution is slightly skewed, the median value is f2 TO = 0.27. 
A 2cr clipping rejects two measurements and gives fi m = 
0.28 ± 0.14. Whereas the mean value turns out to be quite 
stable, it appears to be more difficult to get a stable estimate 
of the error which includes cosmic variance. A still larger 
number of simulations is necessary to improve the accuracy. 
For the comparison with the internal errors given by the KL 
method (see Sect.^) we use the latter more stable estimate, 
keeping in mind that the error could be underestimated by 
about 30%. 

For the linear matter normalization the following formal 
means and standard deviati ons are obtained: as = 0.87 ± 
0.32 (high-peak biasing, Fig. Al upper right), as = 1.20 ± 

"~| lower left), 0.82 ± 0.43 
lT lower right). Note that 



0.66 (Mo & White biasing, Fig. 
(Sheth & Tormen biasing, Fig.L 
for the computation of the means and standard deviations 
of the latter two biasing models only the value s located in 
the high biasing regime are used (see Sect. 3.2.2). The errors 
include cosmic variance. 
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